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Abstract 

Background: Patellofemoral joint replacement is a successful treatment option for isolated patellofemoral 
osteoarthritis. However, results of later conversion to total knee replacement may be compromised by 
periprosthetic bone loss. Previous clinical studies have demonstrated a decrease in distal femoral bone mineral 
density after patellofemoral joint replacement. It is unclear whether this is due to periprosthetic stress shielding. The 
main objective of the current study was to evaluate the stress shielding effect of prosthetic replacement with 2 
different patellofemoral prosthetic designs and with a total knee prosthesis. 

Methods: We developed a finite element model of an intact patellofemoral joint, and finite element models of 
patellofemoral joint replacement with a Journey PFJ prosthesis, a Richards II prosthesis, and a Genesis II total knee 
prosthesis. For each of these 4 finite element models, the average Von Mises stress in 2 clinically relevant regions of 
interest were evaluated during a simulated squatting movement until 120 degrees of flexion. 

Results: During deep knee flexion, in the anterior region of interest, the average Von Mises stress with the Journey 
PFJ design was comparable to the physiological knee, while reduced by almost 25% for both the Richards II design 
and the Genesis II total knee joint replacement design. The average Von Mises stress in the supracondylar region of 
interest was similar for both patellofemoral prosthetic designs and the physiological model, with slightly lower 
stress for the Genesis II design. 

Conclusions: Patellofemoral joint replacement results in periprosthetic stress-shielding, although to a smaller 
degree than in total knee replacement. Specific patellofemoral prosthetic design properties may result in differences 
in femoral stress shielding. 
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Background 

Patellofemoral joint replacement is a successful treatment 
option for isolated patellofemoral osteoarthritis [1,2]. Only 
the patellofemoral joint is replaced, and the femorotibial 
compartments with cruciate ligaments and menisci are 
spared, which probably allows preservation of physiological 
femorotibial joint mechanics. The long-term outcomes of 
patellofemoral joint replacement are related to progression 
of femorotibial osteoarthritis and the need for conversion 
to total knee replacement [3]. Loss of distal femoral bone 
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may compromise the results of such a conversion, and 
therefore needs to be prevented as much as possible. After 
total knee replacement, loss of bone occurs due to the 
stress shielding effect of the femoral component [4,5]. Al- 
though the femoral component of a patellofemoral pros- 
thesis is smaller than in total knee replacement, it is 
unknown whether mechanically induced periprosthetic 
bone remodeling occurs following patellofemoral joint 
replacement. 

Measurements of the periprosthetic bone mineral 
density (BMD) using dual-energy x-ray absorptiometry 
(DXA) demonstrated a 15% decrease in BMD behind the 
anterior flange of the femoral component during the first 
year after Richards II patellofemoral joint replacement 
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[6], but it is not known whether this is due to stress 
shielding. Finite element analyses have been used exten- 
sively in the evaluation of prosthetic load, stress distribu- 
tion in the bone and bone remodeling after total hip and 
knee replacement [7,8]. Some investigators have also 
used numerical models to calculate the stress distribu- 
tion within the patellar components after patellofemoral 
joint replacement [9,10]. However, none of these models 
analyzed a loaded distal femur during a squat to investi- 
gate the effect of the femoral component on the stress 
distribution in the periprosthetic bone. 

Our hypothesis was that patellofemoral joint replace- 
ment results in stress -shielding in the distal femur, but 
to a lesser extent than following total knee replacement. 
The objective of the current study was thus to investi- 
gate the effect of patellofemoral replacement on the 
expected stress distribution in the distal femur eventu- 
ally leading to changes in bone density. For this purpose, 
the patellofemoral joint was modeled in a dynamic finite 
element knee model with and without a patellofemoral 
joint replacement. Furthermore, to investigate the effects 
of different patellofemoral joint replacement designs, we 
compared the Von Mises stresses between the Journey 
PFJ patellofemoral prosthesis, the Richards II patellofe- 
moral prosthesis, and the Genesis II total knee pros- 
thesis in two clinically relevant regions of interest. 

Methods 

Physiological knee joint geometry 

Magnetic Resonance (MR) imaging data (slice thickness 
1.5 mm, pixel dimension 0.43 mm) from an intact 
human right knee cadaver specimen without known 
osteoarthritic changes were manually segmented using 
MIMICS 13 and 3-matic 4.2 software (Materialise, 
Belgium) and reconstructed into a three-dimensional 
model of the osseous and cartilaginous geometries of 
distal femur and patella with patellar tendon and inser- 
tion of the quadriceps tendon (Figure 1). 

Knee joint material properties 

The trabecular bone of the distal femur was modeled as a 
homogenous isotropic linear elastic material (£ = 300 MPa, 
Poisson ratio v = 0.30, density =1 g/cm 3 ) [11,12], while the 
cortical bone was modeled as an orthotropic 2 mm thick 
layer (£1 = 17900 MPa, £2 = 18800 MPa, £3 = 22800 MPa, 
G 23 = 7110 MPa, G 3 i = 6580 MPa, G 12 = 5710 MPa) [11- 
13]. The patella was modeled as a homogenous isotropic 
linear elastic material (£=15000 MPa, Poisson ratio 
v = 0.30, density = 2 g/cm 3 ), with a 2.5 mm thick layer of 
homogenous isotropic linear elastic cartilage (£ = 5 MPa, 
Poisson ratio v = 0.46, density = 1 g/cm 3 ) in contact with 
the femoral trochlea [12,14-16]. The friction coefficient 
between patella and trochlear groove was set at 0.001 
based on experimental data [17]. The patellar tendon 




Figure 1 Physiological model geometry used in this study. 

I J 



was modeled as an isotropic and hyperelastic material 
[16,18], with rubber-like mechanical behavior. 

Knee joint load and constraints 

The simulated motion consisted of a 10s loaded full 
squat (one cycle), starting from 0° until a maximum 
flexion angle of 120°. These settings match the experi- 
mental kinematics simulations performed in a previous 
in vitro analysis on physiological cadaver legs [19-21]. 
The patella model was constrained by fixing the distal 
part of the patellar ligament and applying a quadriceps 
force distributed on the quadriceps insertion on the 
proximal surface of the patella, resulting in the patella 
moving along the trochlear surface of the femur [22]. 
The magnitude and direction of the quadriceps force as 
well as the three-dimensional kinematics of the femur 
were derived from the above mentioned tests on healthy, 
full leg cadaver specimens [21] and were applied using 
the Grood-Suntay coordinate system [23], using as origin 
the midpoint between the condylar center [19-21]. 
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Physiological knee joint finite element model 

The three-dimensional models of femur, patella and pa- 
tellar tendon were imported in commercially available fi- 
nite element analysis software (Abaqus 6.8EF-1, Simulia, 
Dassault Systemes Paris, France). The models were me- 
shed with 2 mm 6-noded triangle elements for the cor- 
tical and cartilage layers, and 2 mm 10-noded 
tetrahedral elements for the trabecular bone, patellar 
bone and patellar ligament. The number of elements for 
each component of the physiological knee model is given 
in Table 1. Two regions of interest (ROI) were defined 
in the femoral bone: an anterior and a proximal ROI. 
The location of the ROIs was defined to fit the same 
regions as used in a previous BMD analysis following 
patellofemoral joint replacement, thus allowing compari- 
son with clinical data (Figure 2) [6]. The ROIs were 
1 cm high in femoral proximo-distal direction and 1 cm 
long in the anteroposterior direction. They spanned the 
entire medio-lateral width. During the dynamic simula- 
tion the average Von Mises stress in each ROI was 
calculated. 

Patellofemoral joint replacement finite element models 

Three different joint replacement designs were considered 
in this study: Journey PFJ patellofemoral joint replacement 
(Smith&Nephew, Memphis, TN, USA), Richards II patello- 
femoral joint replacement (Smith&Nephew, Memphis, TN, 
USA), and Genesis II total knee replacement (Smith& 
Nephew, Memphis, TN, USA). 

The geometries of the Journey PFJ and Richards II pros- 
thetic components were taken from CAD files provided by 
the manufacturer (Figures 3 and 4). The Journey PFJ fem- 
oral component (size medium) and 032 mm patellar com- 
ponent were incorporated in the knee model following the 
manufacturers' instructions, and the Richards II femoral 
and patellar components (size small) were incorporated in 
the model using the previously described surgical technique 
[3]. The Journey PFJ femoral component (oxidized zirco- 
nium or Oxinium™) was modeled as an isotropic linear 
elastic material (£ = 97905 MPa, Poisson ratio v = 0.3, 
density = 6.62 g/cm 3 ) [24], while the Richards II femoral 



Table 1 Number of elements for each component for the 
different finite element models 



Finite element 
model 


Femur 


Patella 


Patellar 
ligament 


Femoral 
component 


Patellar 
component 


Physiological 
model 


326737 


93830 


6914 


/ 


/ 


Journey PFJ 
model 


496078 


26791 


6914 


18396 


13346 


Richards II 
model 


318408 


9375 


6914 


19522 


4636 


Genesis II 
model 


160486 


26791 


6914 


65503 


13346 




Figure 2 Location of anterior (A), and proximal (P) regions of 
interest (ROI) on a lateral radiograph of a right knee following 

Richards II patellofemoral joint replacement. 

k J 



component (CoCr) was modeled as an isotropic linear elas- 
tic material (£=240000 MPa, Poisson ratio v = 03) [25]. 
The UHMWPE patellar component was modeled as a non- 
linear elasto-plastic material (£ = 684.65 MPa, Poisson ratio 
v = 0.45) [19,26,27]. A 2 mm cement layer with linear elastic 
material properties (£ = 3000 MPa, Poisson ratio v = 0.3) 
was modeled between the prosthetic components and the 
cut bone surfaces [28,29]. The frictional coefficient between 
the UHMWPE patellar component and the Journey PFJ 
Oxinium femoral component was set at 0.04 based on ex- 
perimental data [30], and we used a friction coefficient of 
0.08 in the contact zone between the polyethylene patellar 
component and the Richards II CoCr femoral trochlear 
component [30]. The models were meshed with 2 mm 6- 
noded triangle elements for the cortical bone and cement 
layers, and 2 mm 10-noded tetrahedral elements for the 
cancellous bone, patellar bone, femoral component, patellar 
component and patellar ligament. The number of elements 
for each component of the patellofemoral joint replacement 
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models is given in Table 1. We defined an anterior and a 
proximal ROI in the same positions and with the same 
dimensions as in the physiological model During the dy- 
namic simulation the average Von Mises stress in each ROI 
was calculated. 

Total knee replacement finite element model 

The geometries of the Genesis II total knee prosthetic com- 
ponents were taken from CAD files provided by the manu- 
facturer (Figure 5). The Genesis II posterior stabilized total 
knee femoral component (size 5) and 032 mm patellar 
component were incorporated in the knee model following 
the manufacturers' instructions. The Genesis II total knee 
prosthesis was chosen for the total knee model because the 
geometry of the articular surface of the femoral component 
of this prosthesis is exactly the same as the geometry of the 
articular surface of the Journey PFJ femoral component. 
The material properties of the femoral component (Oxi- 
nium™), patellar component and cement were the same as 
in the Journey PFJ model. Also, the contact, friction and 
mesh properties were the same as in the Journey PFJ 
model. The number of elements for each component of the 
Genesis II model is given in Table 1. We defined an anter- 
ior and a proximal ROI in the same positions and with the 
same dimensions as in the physiological model. During the 
dynamic simulation the average Von Mises stress in each 
ROI was calculated. 



Convergence analysis and validation 

A convergence analysis was performed to evaluate the ef- 
fect of mesh quality and density changes on the calculated 
stress distributions. We created 5 mesh files of different 
uniform element sizes ranging from 4 to 2 mm. For each 
mesh size, the contact force, contact area, and contact pres- 
sure were calculated as a function of flexion angle. With a 
sequence of finer meshes, i.e. increasing the mesh density, 
the curves converged to the same solutions. Similarly, we 
ensured that models with different numbers of elements 
converged. 

We compared the numerically determined contact 
area with the contact area obtained using contact pres- 
sure sensors in a knee kinematics simulator experiment 
using the same specimen, with the same load and 
boundary conditions as those used for the numerical 
simulation [31]. The results confirmed the validity of the 
model as shown in Figure 6. Moreover, the values found 
for patellofemoral peak force were in agreement with 
values reported in the literature [22,32,33] . 

Results 

Figure 7 shows the predicted average Von Mises stress in 
the anterior and proximal ROI for the 4 models at discrete 
flexion angles of simulated squat motion. Overall, the aver- 
age Von Mises stress in both ROIs increased with the 
flexion angle. Maximum stresses during squat were reached 
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at 90° flexion angle (2.8-3.8 MPa for the anterior ROI, and 
1.4-1.6 MPa for the proximal ROI). 

Mean stresses in the proximal ROI were similar for 
both patellofemoral joint replacement designs and the 
physiological model, and slightly lower for the Genesis 
II total knee prosthetic design. Between 80° and 120° 
flexion, anterior ROI bone stresses for the Journey PFJ 
design were comparable to the physiological knee, while 
reduced by almost 25% for both the Richards II patello- 
femoral joint and Genesis II total knee joint replace- 
ment designs. 

Discussion 

The aim of the current study was to evaluate the effect of 
patellofemoral joint replacement on the expected stress - 
distribution in the distal femur leading to bone density 
changes. We modeled a patellofemoral joint in a finite 
element knee model with and without a patellofemoral 



joint replacement. In order to investigate the effects of 
different patellofemoral replacement designs, we also 
compared the Von Mises stress distribution induced by the 
patellofemoral load between the Journey PFJ patellofemoral 
prosthesis, the Richards II patellofemoral prosthesis, and 
the Genesis II total knee prosthesis. 

In summary, our dynamic finite element modeling 
demonstrated comparable Von Mises stress in the anter- 
ior region of interest for the Journey PFJ design and the 
physiological knee. The Richards II design demonstrated 
lower average Von Mises stress in the anterior region of 
interest compared to the physiological knee, and the 
Genesis II total knee design demonstrated the lowest 
average Von Mises stress. The average Von Mises stress 
in the proximal region of interest were similar for both 
patellofemoral designs and the physiological model, with 
slightly lower average Von Mises stress for the Genesis 
II total knee prosthesis. 
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Figure 6 Comparison of contact area versus flexion angle for numerical model and knee kinematics simulator experiment. 
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Figure 7 Average Von Mises stress in the anterior and proximal ROI in the 4 models for discrete flexion angles of simulated squat 
motion. 



Stress shielding behind the anterior flange of a patello- 
femoral prosthesis may result in mechanically induced 
bone remodeling with resulting decrease in BMD. Our 
Journey PFJ model predicted no significant stress shield- 
ing, implying that the physiological strains are main- 
tained. However, the Richards II model predicted a 
reduction in average Von Mises stress compared to the 
physiological model. This is in agreement with the re- 
cently reported results of clinical dual-energy X-ray ab- 
sorptiometry (DXA) measurements obtained in 14 
patients [6]. In this prospective 1-year DXA study a 15% 
decrease in BMD was found behind the anterior flange 
of the Richards type II patellofemoral prosthesis. 

The observed differences in the stress shielding effect 
between the Richards II and Journey PFJ prostheses may 
result from differences in both material and geometrical 
properties. The Journey PFJ prosthesis is considered a 
third generation patellofemoral prosthesis. The implant 
is characterized by a broad anatomical trochlear compo- 
nent made of oxidized zirconium (Oxinium™) which has 
a coefficient of friction that is half that of CoCr. There 
are 4 small fixation pegs on the posterior aspect of the 



femoral component. In contrast, the Richards II pros- 
thesis is a first-generation CoCr prosthesis. The non- 
anatomic trochlear component is highly constrained 
with a deep central groove, and the polyethylene patellar 
component has a longitudinal ridge. In addition, the 
femoral component has an approximately 2 cm long 
central fixation peg for added stability. This long central 
fixation peg may in part explain the pronounced stress 
shielding effect observed at 90 degrees of flexion, since 
the patella contacts the femoral component below the 
position of the peg. Furthermore, the Journey PFJ fem- 
oral component is shorter in the sagittal plane than the 
Richards II component. As a result, in deep flexion the 
patellar button contacts the native femoral surface beyond 
90 degrees of flexion while in the Richards II design there 
is still contact between the femoral component and the 
patellar button. 

After total knee replacement, the BMD of the distal 
femur decreases by 16-36% within one year because of 
the femoral components stress-shielding effect [34-40]. 
Finite element analyses have shown that with a bonded 
femoral component, the predicted long-term bone loss 
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would occur at the most distal part of the femur and be- 
hind the anterior part of the prosthesis [4,5]. This is in 
agreement with our findings. 

Only one previous study reported results from a finite 
element knee model in patellofemoral joint replacement 
[10]. Using the geometries of the Low Contact Stress (LCS), 
Vanguard, and Scout patellofemoral prostheses, a three- 
dimensional finite element model of each particular design 
was created. Stress distributions within the patellar compo- 
nents were calculated for three common daily activities: 
walking, ascending stairs, and rising from a chair. Static 
loading was applied with 420 N patellofemoral joint reac- 
tion force in 15° of flexion (walking gait), 1760 N in 45° of 
flexion (stair ascent), and 1950 N in 90° of flexion (chair 
rise). Based on these models, the authors concluded that 
different patellofemoral implant design geometries influ- 
ence the polymer stress within the patellar component. To 
minimize contact and delamination stresses, they further 
stated that contemporary designs should employ a broad 
trochlear groove to maximize congruent patellar compo- 
nent contact. No attempt was, however, made to evaluate 
for altered stress distributions in the periprosthetic distal 
femur. In an experimental setup using tri-axial strain 
gauges in synthetic femurs before and after Journey PFJ 
patellofemoral joint replacement, Meireles et al. determined 
the strain shielding effect in the distal femur [41]. A strain 
shielding effect was found with static loading of the patello- 
femoral joint in 12, 50, and 90 degrees of flexion. The 
results showed a reduction in strain in the medial and distal 
regions of the femur when deep bending occurred, with 
higher values of strain in the anterior region proximal to 
the prosthesis. This distal diaphyseal increase in strain has 
also been noted in finite element models with total knee re- 
placement [5]. 

The present study is the first to specifically evaluate for 
the stress-shielding effect of a patellofemoral prosthesis 
using finite element modeling. However, correct interpret- 
ation of results obtained from numerical models requires 
careful consideration of several important issues [42]. We 
employed a numerical model that has been both verified 
and validated using the same procedure as described by 
Catani et al. [31]. Moreover, the biological findings based 
on clinical DXA measurements behind the anterior flange 
of the Richards II patellofemoral prosthesis [6] match the 
findings of this finite element analysis. Of course, validation 
will only be complete by adding clinical DXA measure- 
ments in the same ROIs for the other types of implant that 
were considered in this simulation. This information is not 
available yet but it is the subject of an ongoing study. 

Finite element modeling is, however, subject to limita- 
tions due to inherent uncertainties concerning geometry, 
load situation, and material properties. We performed our 
modeling using a dynamic loaded squat instead of a gait 
cycle, which is a more frequent activity. Squatting is a 



motor task which produces higher joint contact forces and 
therefore higher stresses and strains in the bone compared 
to gait [43,44]. Moreover, as squatting gives rise to higher 
stresses, it should also induce a larger difference in stress 
between the physiological and replaced knee. Finally, the 
analysis of a squat, up to 120°, was preferred since it reaches 
a larger range of motion compared to gait which is usually 
limited to 60° of flexion [45]. Nevertheless, analysis of other 
activities will probably lead to similar relative results but in 
a smaller range of knee flexion and with smaller differences 
in periprosthetic Von Mises stress values. Although we 
employed a dynamic simulation, we did not consider visco- 
elastic properties of the materials. Because the knee flexion 
from 0° to 120° was performed in 10s, the time dependent 
effect was assumed negligible. Furthermore, in the simu- 
lated activity, the patellar tendon was always loaded during 
the entire cycle and no unloaded situation or hysteresis was 
present nor simulated. 

A further limitation of our study was that we evaluated 
the stress -shielding effect of three different joint replace- 
ment designs from one manufacturer (Smith&Nephew, 
Memphis, TN, USA). Prosthetic designs from other manu- 
facturers differ with respect to material and design proper- 
ties. Inclusion of other designs in our FE modeling would 
probably have shown differences in the amount of stress 
shielding, and not including these designs therefore limits 
the scope of our study. 

Finally, we did not consider potential failure of the bone- 
implant interface due to elevated stress levels. Although less 
stress shielding means better load transfer to the bone, 
which is desirable, it may also imply higher loads in the 
interface region and an increased risk for failure of the fix- 
ation that was not considered in the present analysis. How- 
ever, this limitation is not clinically relevant as loosening of 
cemented femoral components is not an issue with isolated 
patellofemoral joint replacement [3]. 

Conclusions 

Dynamic finite element analyses of knee models with a 
patellofemoral joint replacement predict a decrease in aver- 
age bone stress behind the anterior flange of the femoral 
component. This reduction was more pronounced in the 
Richards II design than in the Journey PFJ design, and may 
be related to specific design properties. 
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